Dynamical stabilization of solitons in cubic-quintic nonlinear Schrodinger model 
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We consider the existence of a dynamically stable soliton in the one- dimensional cubic-quintic 
nonlinear Schrodinger model with strong cubic nonlinearity management for periodic and random 
modulations. We show that the predictions of the averaged cubic-quintic NLS equation and modified 
variational approach for the arrest of collapse coincide. The analytical results are confirmed by 
numerical simulations of one-dimensional cubic-quintic NLS equation with rapidly and strongly 
varying cubic nonlinearity coefficient. 
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Collapse phenomena are observed in many areas of 
physics: self-focusing of intense laser beams, Langmuir 
waves in plasma, collapse of the Bose-Einstein conden- 
sates (BECs) with attractive interactions etc. The non- 
linear Schrodinger equation (NLSE) with cubic nonlin- 
earity used to describe these systems has stable solu- 
tions in the one-dimensional (ID) case, when the dis- 
persion and nonlinearity effects can effectively balance 
each other. In two and three dimensions the focusing 
nonlinearity overcomes the dispersion and the blow-up 
phenomenon occurs 0]. 

Few mechanisms for the arrest of collapse have been 
suggested. Among them can be mentioned the dis- 
persion and nonlinearity management methods 
120- 0- E3- 13 The analysis based on the variational ap- 
proach, method of moments and numerical simulations 
showed that the nonlinearity management method is ef- 
fective to suppress collapse in the scalar and vector 2D 
NLSE with focusing cubic nonlinearity. For the 3D cu- 
bic NLSE with nonlinearity management theoretical pre- 
dictions and numerical results do not bring a clear and 
definitive picture 0, fl(i| , so more analytical and numeri- 
cal work is necessary. 

In this Rapid Communication we investigate the phe- 
nomenon of arrest of collapse by using the strong cubic 
nonlinearity management scheme in the ID cubic-quintic 
(CQ) NLSE. The CQ NLSE with nonlinearity manage- 
ment presents practical interest since it appears in many 
branches of physics such as nonlinear optics and BEC. 
In nonlinear optics it describes the propagation of pulses 
in double-doped optical fibers , in BEC it models the 
condensate with two and three body interactions 
In optical fibers periodic variation of the nonlinearity 
can be achieved by varying the type of dopants along 
the fiber. In BEC the variation of the atomic scattering 
length by the Feshbach resonance technique leads to the 
oscillations of the mean field cubic nonlinearity 



The CQ NLSE when the cubic term is equal to zero 
is the critical quintic NLSE. The quintic Towncs soli- 
ton is an unstable solution of the quintic NLSE 1]. In 
this work we consider the configuration with the rapid 
and strong periodic modulation in time of cubic nonlin- 
ear interaction. This type of modulations corresponds to 
the management applied to the nonlinearity with a lower 
power and has never been studied. We first apply a vari- 
ational approach to the averaged NLSE. The averaged 
equation for the ID cubic NLSE in the case of st rong; 
nonlinearity management has been derived in [Til Il6| . 
but the presence of the quintic term dramatically changes 
the picture as this term would normally lead to collapse. 
We also propose a modified variational approach for the 
managed CQ NLSE where the ansatz is designed to take 
into account the fast self-phase modulation due to the 
cubic nonlinearity management. These two approaches 
predict the stabilization of the quintic Townes soliton. 

We consider the CQ NLSE 



iu t + u xx + j(t)\u\ 2 u + x\u\ 4 u = . 



(1) 



with an attractive quintic nonlinearity x > 0. The time- 
varying cubic coefficient j(t) possesses an average value 
70 and a fast varying part 71 

7(0=7o + -7i(-), (2) 

e e 

where e <C 1 corresponding to strong and rapid man- 
agement. Here 71 can be either a periodic function 
or a stationary random function. We first address the 
case of a periodic management 71 (t + 1) = 71 (r) and 
Jp 1 71 (t)oIt — 0. We introduce Ti(t) = JJ" 71 (s)ds which 
is also a zero-me an p eriodic function. Following the same 
procedure as in [lfj, we can average the CQ NLSE over 
fast variations and show that the solution takes the form 



u(t, x) 



w(t, x) + EW\(t, -,x) + 

e 



x exp 



e 
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where w is solution of the averaged CQ NLSE 
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iw t + w xx + 7o|w| 2 w + x\w\ Aw 
+<J 2 {l(\w\ 2 ) x } 2 + 2\w\ 2 (\w\ 2 ) xx }w = 0, (3) 

and a 2 = f^T^sfds. The averaged CQ NLSE has a 
Hamiltonian form, with the Hamiltonian 



H 



(\Wx\ 



70 | ,4 

2>l 



3 1 1 



<t 2 [(M 2 )J 2 H 2 )^. 



(4) 

We first prove that the solution of Eq. (J3J) cannot collapse 
because its supremum norm can be a priori bounded. 
This proof is essentially based on the Sobolev inequal- 
ity 11/11*. < C||/|| 2 ||/ x || 2 where we denote ||/|| p = 

\J \f(x)\ p dx~\ 1 P for p G (l,oo) and ||/||oo is the (essen- 
tial) supremum of |/|. We first apply this inequality with 
f = w: 

Nil < C\\ W \\l\\w\\l < CHIalKlla < C\\w x \\ 2 , (5) 

where C stands for a constant that may change from line 
to line and we have used the fact that ||w||2 is constant. 
Next we apply the Sobolev inequality with / = v := \w\ . 

IMIS < NltllH!2 = IIHI 4 oo /3 IIHl2<^IIHl2 /3 IKIl2 /3 



< C8\\v\ 



4/3 



,4/3 



The last estimate holds true for any 8 > 0, so by choosing 
s = \\v\\t /3 /(2C) and noting that ||t>||| = ||io||g, we get 

\Mt<C\\v x \\ 2 . (6) 

Substituting (l- r >l(i[) into J2J| and using the fact that 
/[(IH 2 )a:] 2 |H 2 ^ x — (4/9)||ua;||2, we can finally write 

H>\\w x \\%- l0 C\\w x \\ 2 + (A/9)a 2 \\v x \\ 2 - X C\\v x \\ 2 , 

which shows that ||u> x || and \\v x \\, and thus ||ttf||oo are 
uniformly bounded, which prevents the solution of the 
averaged CQ NLSE from collapsing. However, this re- 
sult does not show that the solution of Eq. © does not 
collapse for fixed e. 

The next step consists in applyin g th e variational ap- 
proach to the averaged CQ NLSE [H E3- The vari- 
ational ansatz for the solution is chosen as the chirped 
function 



w{t,x) = A(t)Q 



i(t) 



exp [ib(t)x 2 + i<f>(tj\ , (7) 



with a given shape Q. Following the standard procedure, 
we substitute the ansatz into the Lagrangian density gen- 
erating Eq. 10 and calculate the effective Lagrangian 
density in terms of A, a, b, <j> and their time derivatives. 
The evolution equations for the parameters of the ansatz 
are then derived from the effective Lagrangian by using 
the corresponding Eulcr-Lagrange equations. In particu- 
lar this approach yields a closed-form ordinary differen- 
tial equation (ODE) for the width a: 



where the effective potential is 

6£>i - 2D 3X N 2 D 2l pN 8D 4 a 2 N 2 
U av (a) = — 2 _ + _._, ( 9 ) 

the total mass is N = J \u\ 2 dx and the effective parame- 
ters are 



Di = 
D 2 = 
D 3 = 

Da = 



jQ'(x) 2 dx 
J x 2 Q 2 (x)dx ' 

jQ(x) A dx 
(jQ 2 (x)dx) (Jx 2 Q 2 (x)dx) ' 

jQ(x) 6 dx 
(jQ 2 (x)dx) 2 (Jx 2 Q 2 (x)dx) ' 

jQ(x) 4 Q'(x) 2 dx 
(jQ 2 (x)dx) 2 (Jx 2 Q 2 (x)dx) ' 



Another approach is possible for the CQ NLSE with 
cubic nonlinearity management Q). It consists in apply- 
ing a specific variational approach designed to capture 
the fast self-phase modulation induced by the fast non- 
linearity management. The variational ansatz is sought 
in the form 



'(('••<•) = Mi)() ( -iL ) exp \ib(t)x 2 + i<p(t)] 



x exp 



l T 1 ( t -)A 2 (t)Q 2 (^-) 
e a(t) 



(10) 



Substituting this ansatz into the Lagrangian density gen- 
erating Eq. (JTJ we get the system of ODEs 



a, 
hi 



Aab 



3?r l( i), 



a 4 
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D2I0N 

4a 3 



D3XN 2 
3a 4 



8D 4 N< 



(11) 

LhNb t 

(12) 
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Next we perform an averaging of this system of ODEs, 
which yields exactly Eq. (|SJ • This analysis shows that the 
two different approaches yield the same effective system 
for the soliton parameters, which strengthens its validity. 

We now focus our attention to the particular case 70 = 
0. We first choose the shape function Q corresponding 
to the quintic Townes soliton Q(x) = l/^/cosh(a;), which 
gives the values D x = l/(2vr 2 ), D 2 = 8/tt 4 , D 3 = 2/tt 4 , 
and D4 = l/(87r 4 ). In absence of cubic nonlinearity man- 
agement a = 0, there exists an infinity of fixed points a c 
(for which U' av {a c ) — 0) if the total mass is equal to the 
critical mass N c = 7T\/3/(2y / x). As it can be checked 
from (JSJl the Townes soliton is not stable in the sense 
that it collapses if N > N c and it spreads out and van- 
ishes if N < N c , so this theoretical solution of the quintic 
NLSE with N = N c cannot be observed in practice. In 
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FIG. 1: Rms pulse width versus time. The initial pulse has 
a Townes shape with mass N = 1.02AT C ~ 2.77 and radius 
ao = a c — 0.88 (corresponding to an initial rms pulse width 
a-rmsfi = 1.38, solid line) or ao = 1.0 (a rmSj o = 1.57, dashed 
hnes). Picture a: collapse at finite time in absence of nonlin- 
earity management. Picture b: a periodic cubic nonlinearity 
management j(t) = 10cos(50t) is applied and the pulse rms 
width is plotted till time t — 350. 



presence of strong cubic management a > 0, we get the 
existence of a unique fixed point a c if N > N c , with 



1 



(13) 



The linear stability analysis shows that this fixed point 
is stable and that the oscillation period is 



2a 2 N 2 



[(N/N c ) 2 - l]3/2 



(14) 



We are especially interested in solutions whose masses 
are just above the critical mass, since our main goal is to 
prove that the Townes soliton can be stabilized by cubic 
nonlinearity management. In these conditions, the stable 
soliton width is rather large ~ 1/(N — Nc) 1 / 2 , and the 
soliton oscillation period is very long ~ l/(N — N c ) 3 / 2 . 

In the following numerical experiments, we take x = 1, 
so the critical mass is N c ~ 2.72, and we apply the man- 
agement 7(f) = 10cos(50i), which gives a = l/(5y2)- 
We solve the CQ NLSE by a pseudo-spectral method 
starting from a Townes soliton-shape with a mass N > 
N c and radius a . For N = 1.02N C (Fig. [TJ, the theo- 
retical fixed point and oscillation period are a c = 0.88 
and T c = 38. If we choose oo = a c for the input pulse 
width, then we observe in the numerical simulations that 
the mean-square radius of the pulse is almost constant, 
which shows that this solution is stable. If we choose 
ao = 1 > a c , then we observe slow oscillations around a c 
with a period very close to T c . These two observations 
demonstrate that the values a c and T c for the soliton pa- 
rameters predicted by the variational analysis are very 
accurate. 

We have repeated the same experiments with different 
initial masses, to check the validity of formulas H13I14JI . 
For N = 1MN C (resp. N = 1.08N C ), the theoretical 
fixed points and oscillation periods are a c = 0.63 and 
T c = 14 (resp. a c = 0.42 and T c — 5.1). As it can be seen 



FIG. 2: Rms pulse width. Picture a: N = 1.0AN C , a = 0.63 
(a T ms,o = 0.99, solid line), ao = 0.75 (a rms ,o = 1.18, dashed 
lines). Picture b: N = 1.08iV c , a = 0.5 (a rms , = 0.79, solid 
line). 



in Fig. |2)3, when the initial mass is significantly larger 
than the critical mass, the pulse first experiences a strong 
distortion and emits radiation. The resulting pulse is 
still stable, its mass is still well above the critical mass, 
but the fixed point and oscillation period are not well 
predicted by the variational analysis. If the initial mass is 
even larger, then this first phase leads to collapse, which 
shows that the stabilizing effect of the cubic nonlinearity 
management -f(t) = 10cos(50i) is only efficient when the 
initial mass is in the range [N c , l.liVj. 

We now prove the robustness of the solution of the 
CQ NLSE driven by cubic nonlinearity management with 
respect to the initial pulse shape. More exactly, we 
show that a stable solution can be obtained with an ini- 
tial pulse profile which is significantly different from the 
Townes profile. The important conditions that have to be 
satisfied by the initial pulse is that its mass should be just 
above the critical mass corresponding to the input pulse 
profile, and that its initial radius should be chosen in the 
vicinity of the fixed point. These conditions are imposed 
by the analysis of Eq. (JHJ , and they are confirmed by nu- 
merical simulations. Let us consider an arbitrary pulse 
shape Q. The critical mass is then N c — y (3D\) / (D^x) 
and the fixed point in presence of nonlinearity manage- 
ment is a c = 2V2V 'Da/ DIN ~/^(N/N c ) 2 - 1, with the 

period T c = (4V27r J D 4 /^ 1 /2 )ff 2 A r2 /[(A r /iV c ) 2 - 1] 3/2 - A 
stable solution of the managed CQ NLS equation can be 
obtained by injecting a pulse with shape Q not too far 
from the Townes soliton, with a mass just above A^ c and 
a radius close to a c . In the case of the Gaussian ansatz 
Q{x) — exp(— x 2 /2), we have D\ = 1, D 2 = 
D 3 = 2/(v / 3tt), and D 4 = 1/(3V3tt). We report in 
Fig. the results of numerical simulations carried out 
with x — 1 so the critical mass is N c ~ 2.86. We use the 
same management as above. In a first phase (i < 20), 
the pulse emits radiation and its shape converges to the 
one of the Townes soliton. A similar phenomenon has 
been observed in numerical simulations carried out in 
[8j for the 2D cubic NLSE with cubic nonlinearity man- 
agement. After this transition period, the soliton width 
experiences oscillations around the fixed point with an 
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FIG. 3: Dynamics of a Gaussian pulse with initial mass N = 
1.02jV e = 2.92 and radius ao = a c = 1.4 (a rms ,o = 0.99). 
Picture a: pulse profile at time t = 350 (solid line) compared 
to the input Gaussian profile (dashed lines). Picture b: rms 
pulse width versus time. 
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FIG. 4: Rms pulse width versus time. The initial pulse has 
a Townes shape with mass TV = 1.02iV c ~ 2.77 and radius 
ao = a c = 0.88 (a rms ,o = 1.38, solid line) or ao = 1.0 
(arms, o = 1.57, dashed lines). We apply a random cubic non- 
linearity management *y(t) = 10m(50£) where m is zero-mean 
stationary random process with a unit coherence time com- 
pensated by a pinning scheme. 

oscillation period close to the theoretical value. 

We finally study the stabilization induced by random 
nonlinearity management. Accordingly we now consider 
that 71 is a zero-mean stationary random process. Con- 
sidering Eqs. i|llll2fl . it can be seen that the important 
process is actually Ti(t) = J Q ji(s)ds, and it is critical 
that this process does not grow in a diffusive manner, 
which would mean that cubic nonlinearity accumulates. 
This condition is fulfilled if 71 is the derivative of a sta- 
tionary random process, or if we apply a pinning scheme. 
This technique was first introduced by j 1 Qj| to compensate 
for accumulated fiber dispersion. The periodic insertion 
of additional pieces of fiber with well controlled lengths 
and dispersion values was found to prevent from pulse 
deterioration. The pinning method can be applied to 
compensate for accumulated cubic nonlinearity as well. 
In Fig. 01 we show that random cubic nonlinearity man- 
agement stabilizes a Townes soliton in a manner similar 
to a periodic management. However, we can detect a very 
slow spreading out, whose origin can be explained by the 
imperfect compensation of the accumulated nonlinearity 
by the pinning scheme. 

In conclusion we have analyzed the stabilizing role of 
the strong management of the cubic nonlinearity in the 



ID cubic-quintic NLSE. We have proved that the aver- 
aged CQ NLSE, in a dramatic distinction from the non- 
modulated system, supports stable solutions beyond the 
critical mass. In particular the quintic Townes soliton, 
which is unstable in the non-modulated system, becomes 
stable in presence of strong nonlinearity management if 
its width lies in the vicinity of some fixed point. This 
fixed point and the associated oscillation frequency are 
well predicted by the variational approach applied to the 
averaged CQ NLSE. We have also developed and applied 
a modified variational approach to the CQ NLSE with 
strong management, which gives the same values for the 
parameters of the stabilized soliton. We have checked 
that a random management is also stabilizing, if we take 
care to select a random process whose cumulative values 
are small. These results can be applied for the search 
of nonlinearity managed solitons in double-doped optical 
fibers and BECs with three-body interactions. 
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